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INTRODUCTION 


Atmospheric  dispersion  models  are  mathematical  expressions  that  combine 
source  emissions  with  meteorological  parameters  to  produce  air  quality 
estimates  at  specified  receptor  locations.  At  airports  where  many 
sources  and  receptors  are  involved,  refined  models  such  as  the  Airport 
Vicinity  Air  Pollution  model  (AVAP)(1)  are  used  to  determine  air  quality. 
However,  where  few  sources  and  receptors  are  involved,  screening  models 
are  very  attractive  for  identiving  the  need  for  further  analysis  with 
the  more  refined  models.  This  report  describes  one  of  these  screening 
models,  SIMPLEX  "A". 

This  report  describes  the  mathematical  basis  for  the  model,  lists  the 
program,  and  explains  the  steps  taken  to  compute  pollution  dosage. 

Special  program  features  are  described  and  two  sample  problems  are 
solved . 


The  experienced  user,  who  is  primarily  concerned  with  running  a  specific 
problem,  may  bypass  the  descriptive  sections  of  this  report  and  proceed 
directly  to  the  "Sample  Problem-Program  Operation"  section  on  page  5. 

MODEL  DESCRIPTION 


SIMPLEX  "A",  which  has  been  programmed  for  the  Hewlett  Packard  67  and  97 
desk  calculators,  addresses  emissions  during  takeoff.  Additional  SIMPLEX 
models  are  being  developed  to  determine  the  air  quality  impact  from 
taxiing  and  queueing  aircraft  as  well  as  from  ground  vehicles  at  the 
airport.  The  model  is  particularly  useful  at  small  airports  and  at 
those  airports  having  only  a  few  dominant  sources. 


SIMPLEX  "A"  uses  the  same  Gaussian  formulation  employed  in  many  of  the 
refined  models  listed  in  the  Environmental  Protection  Agency’s  (EPA) 
guidelines  on  air  quality  models.  It  accomplishes  its  function  by 
simplifying  many  of  the  detailed  features  of  the  more  refined  models. 

The  model  is  an  integrated  puff  model  for  an  accelerating  point  source. 
Downwind  receptors  are  assumed  to  be  at  ground  level  (z=0)  and  receive 
pollution  doses  from  each  emission  puff.  Figure  1  describes  the  source- 
receptor  geometry  where  the  dose  from  each  emission  puff  is  summed  at  a 
receptor  to  give  a  total  dose  due  to  a  complete  takeoff  event.  Concentrations 
are  measured  in  parts  per  million  (ppm)  of  Nitrogen  Dioxide  (NO^)  where 
the  complete  conversion  of  Nitrogen  Oxides  (NO  )  to  NO2  is  assumed.  In 
cases  where  Carbon  Monoxide  (CO)  concentrations  are  required,  NO^ 
calculations  can  be  factored  appropriately.  The  total  dose  at  point 
x,y,0  is  given  by  the  equation: 


Q 


tt 


SYMBOL 


—  CaP 
"'TV  U 
DEFINITION 


F‘(t.  )' 

L  rr 


L  ^ '  *  -  * 


I  T 

UNITS 


J 


(2) 


V  =  receptor  exposure  of  dose 

*  =  downwind  distance  in  the  direction 

of  the  mean  wind 

Y  =  crosswind  distance 


ppm-sec.  (NO^) 

meters  (m) 
m 


i 

height  above  ground  level 

m 

standard  deviation  of  plume  concen¬ 

*T 

tration  in  the  vertical  direction 

m 

V 

standard  deviation  of  plume  concen¬ 

tration  in  the  crosswind  direction 

m 

c 

wind  speed 

in/  sea* 

<V 

total  emissions  during  an  emission 

re  1  ease 

grams 

H  = 

effective  heighL  of  Lin  iss  ions 

m 

'1  hi’  program  is  printed  out  in  Figure  2.  The  registers,  labels,  inputs 
and  outputs  aro  listed  in  K  i}',un  s  1  and  . 

SFF.C  1A1.  PKOCIUM  KKATl-RKS 

Standard  Deviations  *>  f  Plume  Concent  ra  t  ion  (sigma  (jT)  ) 


A  subprogram  was  employed  to  tic  term  ine  the  standard  deviation  at  pl:r\. 
concon t ra t ion  in  the  horizontal  (crosswind)  and  vertical  directions. 
This  subprogram  was  based  upon  the  assumption  that  pollution  disperses 
according  to  the  power  law  expression: 

2  (J~  -  Kx^  or,  in  straight  line  form 

3  bog  3T-  b  Log  x  +  Log  K 


The  exponent  "b"  governs  the  rate  of  pollutant  dispersion  and  tile  coot  fir  i 
"K"  depends  upon  atmospheric  stability. 

Analysis  of  the  Pasquil 1 /C if ford  curves*  (  t)  used  in  most  dispersion 
models  shows  that  for  stability  classes  "R"  through  "F",  single  straight 
lines  are  approximated  whenG”v  andCTz  values  are  plotted  logarithmically 
against  downwind  distance  up  to  a  source-receptor  distance  of  1000 
meters.  It  is  also  seen  that  these  straight  1 ines  have  the  same  slope. 


With  the  realization  that  as  a  fund  ion  of  tour  stability  classes  can 
be  represented  as  single  straight  lines  with  the  same  slope  (0.9), 
equation  2  can  be  rewritten  as: 
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CTV,  =  K..,x 
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when  x  does  not  exceed  1000  meters 


The  values  of  K  and  K .  which  are  listed  in  the  program  printout  (Figure 
2)  were  obtained  by  solving  equations  4  and  S  for  and  K,}  after 
substituting  values  for  x,n7"v  a  nd  (T  z. 

The  program  not  only  allows  for  calculation  of  the  (Tv  and  *J~z  values  in 
equations  4  and  *>  but  also  has  provisions  to  inuit  valuer  for  initial 
sigmas  «To  andflf'  )  in  order  to  .u-coimt  for  the  enhanced  dispersion 
caused  by  the  hot,  high  velocity  jet  exhaust.  (This  enhanced  dispersion 
is  discussed  in  (4),  (5)  and  (b)).  The  final  dispersion  parameter  (Ol) 
is  automatically  calculated  in  the  program  bv  summing  and  squaring 
initial  and  ca l cu 1 ated ^ 1 s  and  then  extracting  the  square  root. 


*  The  values  from  the  Pasqu i 1 1 /Gi f ford  curves  were  used  per  se  and  were 
not  adjusted  for  averaging  times  different  from  those  stated  in 
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Because  of  technical  difficulties,  it  lias  not  been  possible  to  determine 
Cfo  values  from  measurements  ta*en  at  airports  during  high  thrust  airplane 
takeoff.  However,  plume  measurements  have  been  made  during  low  thrust 
operations  at  Dulles  (4)  and  Lor,  Angeles  International  (7)  airports. 
Average  values  from  measur ements  taken  at  these  two  airports  (8  meters 
for  /T°Z  and  lb  meters  for  <7o  )  are  !  ncorpora  ted  in  the  model. 

Plume  Height 


Because  of  the  lack  of  experimental  data  to  support  plume  rise  theories 
for  taking-off  aircraft,  special  plume  rise  algorithms  have  not  been 
incorporated  into  the  model.  While  research  is  planned  in  this  area, 
until  this  research  is  completed,  it  was  assumed  that  the  plume  height 
was  at  least  as  high  as  the  airplane  engines.  An  average  value  for  this 
parameter  is  four  meters  for  airplanes  operatng  at  a  typical  large 
airport . 

Stability  Class 

Pasquill /Gif ford  stability  classes  "B",  "D",  and  "E"  are  expected 

to  prevail  at  the  airport  during  the  times  of  air  quality  assessment. 

Turner  (3)  gives  a  detailed  description  of  the  characteristics  of  each 
stability  class.  A  particular  stability  class  is  identified  by  a  range 
of  wind  speed,  solar  radiation  intensity,  and  cloud  cover.  Values  for 
these  parameters  can  be  obtained  from  local  National  Weather  Service  or 
observer  personnel. 

Winds 

The  coordinate  system  is  oriented  to  the  runway  on  which  the  aircraft 
are  assumed  to  be  operating.  Since  aircraft  usually  take  off  into  the 
wind,  wind  angles  are  measured  only  from  0  to  90  degrees  on  either  side 
of  the  runway.  For  example,  a  zero  degree  wind  would  blow  directly  down 
the  runway;  a  90  degree  wind  would  blow  perpendicular  to  the  runwav. 

Accelera  t ion 

Aircraft  performance  manuals  can  be  used  to  determine  acceleration 
during  takeoff.  However,  the  program  has  been  structured  to  accept  an 
average  takeoff  acceleration  and  performance  information  should  be 
adjusted  to  average  acceleration  values. 

Emission  Tail 

During  the  operation  of  a  jet  engine,  the  high  velocity  of  its  exhaust 
gases  creates  an  emission  tail  which  can  extend  for  a  considerable 
distance  behind  the  aircraft  (Figure  5).  This  tail  is  simulated  bv 
assuming  a  value  for  its  length  and  by  assuming  a  finite  number  of 
points  along  the  tail  at  which  emissions  are  considered  to  be  released. 
Observations  of  the  tail  length  of  a  number  of  aircraft  enabled  the 
estimation  of  a  nominal  value  for  this  tail  length  and  parametric  calculation 
of  concentrations  at  receptors  enabled  the  selection  of  a  minimum  number 
of  emission  release  points  that  would  provide  reasonable  model  accuracy 
for  those  wind  angles  that  permit  a  sweep  of  pollution  over  a  receptor 
during  the  takeoff  run. 
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The  program  assumes  a  225  meter  emission  tail  with  three  emission  release 
points  located  75  meters  apart  in  the  tail  (see  Figure  5).  The  model  is 
programmed  to  index  the  emission  starting  point  75  meters  further  down 
the  tail  after  each  iteration  sequence  is  completed.  The  first  and  last 
points  in  the  tail  are  37.5  meters  from  the  ends  of  the  tail.  (The  tail 
starts  at  the  exit  plane  of  the  engine.) 

Vertical  Dispersion  Lid 

Calculations  under  a  variety  of  assessment  conditions  showed  that  a  lid 
on  vertical  dispersion,  i.e.  an  inversion  "cap",  had  an  insigni f irant 
effect  on  concentration  at  the  short  downwind  distances  (less  than  1000 
meters)  employed  in  assessing  aircraft  pollutant  impact.  An  algorithm 
to  account  for  this  phenominon  is,  therefore,  not  included  in  the  program. 

Iteration  Interval 

From  past  program  use,  a  one-second  iteration  time  is  recommended. 

Using  this  iteration  time  interval,  the  dose  calculation  can  he  completed 
in  less  than  15  minutes  at  a  source  receptor  distance  of  300  meters. 

Dosage  Output 

Total  dosage  is  printed  out  at  the  end  of  each  iteration  sequence.  The 
program  is  stopped  when  the  dosage  reaches  a  maximum  value.  Output 
units  are  parts  per  million-seconds  ( ppm-seconds) .  To  determine  the 
average  concentration  over  a  one-hour  time  period  (for  compatabilitv 
with  a  particular  short  term  standard)  the  dosage  must  be  divided  bv 
3600  seconds. 

SAMPLE  PROBLEM 

The  step  by  step  procedure  for  solving  the  sample  problems  is  described 
in  this  section.  While  this  procedure  is  structured  for  a  single  aircraft, 
the  same  procedure  can  be  used  for  anv  number  of  aircraft  bv  treating 
them  as  one  large  aircraft. 

Preparation  of  Data  For  Program  Execution 

From  Figure  1  it  is  seen  that  the  Case  1  receptor  is  located  337.5 
meters  downwind  from  the  aircraft  (as  measured  along  the  runwav)  and  200 
meters  abeam  of  the  runway  centerline.  The  Case  2  receptor  is  located 
262.5  meters  upwind  of  the  aircraft  and  100  meters  abeam  to  it. 

The  objective  of  this  problem  is  to  determine  the  air  qualitv  impact  of 
747  NOx  (‘missions  (reported  as  NO2)  during  takeoff.  During  this  takeoff, 
it  was  assumed  that  a  5-meter  per  second  wind  was  blowing  at  30  degrees 
to  the  runway  centerline  and  that  Pasquill/Gi f ford  stability  class  MEn 
prevailed.  The  747  was  assumed  to  have  a  constant  takeoff  acceleration 
of  1.3  meters  per  second  per  second.  All  imputs  are  listed  in  Figures  6 
and  7. 
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The  following  procedure  was  used  in  solving  the  problems: 

Source  emissions  were  obtained  from  AP-42  supplement  10  (8),  where  747 
NO^  emissions  are  listed  at  215.3  kilograms  per  hour  per  engine  or  60 
grams  per  second  per  engine.  Since  the  747  has  four  engines*  the  total 
emission  rate  was  240  grams  per  second.  To  accommodate  the  three  emission 
release  points  in  the  "tail"  (see  Figure  5)  this  rate  was  divided  by  3 
to  reduce  its  value  to  80  grams  per  second  per  "tail"  release  point. 
Selecting  an  iteration  time  of  one  second  and  multiplying  it  by  the 
emission  rate  results  in  the  release  of  80  grams  of  NO^  per  puff. 

<foz  and  Toy  values  of  8  and  16  meters  respectively  were  selected  from 
the  Standard  Deviation  of  Plume  Concentration  section  of  this  report  and 
a  plume  height  of  4  meters  was  selected  from  the  plume  height  section. 

The  beginning  time  was  set  at  zero  by  inputting  the  iteration  time  (one 
second)  and  assigning  a  negative  sign  to  it.  The  Case  1  receptor  is 
downwind  of  the  aircraft  giving  it  a  positive  sign  (see  Figure  1).  The 
Case  2  receptor  is  upwind  of  the  aircraft  giving  it  a  negative  sign. 

The  airplane  to  receptor  distance  is  converted  to  a  ’’tail"  to  receptor 
distance  (at  the  first  "tail"  emission  point,  see  Figure  5)  by  subtracting 
37.5  meters  from  the  former  to  uniformly  space  the  three  tail  release 
points  over  the  225  meter  tail  length.  The  resulting  distance  between 
the  receptors  and  the  first  point  in  the  emission  tail  is  +300  meters 
feu*  Case  1  and  -300  meters  for  Case  2. 

Prog  ram  Opera  t ion 

Load  the  Program 

Before  loading  the  program  the  ’’on-off”  switch  should  be  in  the  ”on” 
position  and  the  "run-program”  switch  should  be  in  the  "run”  position 
(for  the  HP-97  the  "trace-manual -norm"  switch  should  be  in  the  "manual” 
position).  The  program  can  then  be  loaded  into  the  calculator  by  first 
pushing  the  number  1  end  of  the  magnetic  tape  strip*  into  the  slot  in 
the  upper  left  hand  portion  of  the  HP-97  calculator.  (On  the  HP-67 
calculator  the  slot  is  located  on  the  right  hand  side.)  When  the  strip 
comes  out  the  other  side,  turn  it  around  to  the  number  2  end  and  push  it 
through  the  slot  a  second  time.  The  program  is  now  loaded  into  the 
calculator  and  the  tape  strip  which  has  come  out  the  hack  of  the  calculator 
can  be  stored  in  the  horizontal  slot  just  under  the  calculator  switches. 

input  Data 

Inputs  for  the  Case  l  and  Case  2  problems  are  listed  in  Figures  6  and  7 
and  a  printout  of  the  results  is  listed  in  Figures  8  and  9.  The  Case  1 
problem  is  solved  bv  first  entering  the  values  for  the  six  input  parameters 
1 isted  in  Figure  6  into  the  Primary  Register  by  depressing  the  following 
keys:  80  ST0  0  8  STO  2  4  STO  B  30  STO  D  1  STO  E  5  STO  I .  Any  input 
errors  can  be  erased  by  depressir;  the  CLx  key  or  bv  turning  off  the 
calculator,  restarting  it  and  reloading  the  program. 


*  This  tape  strip  can  be  obtained  bv  contacting  the  Federal  Aviation 
Administration;  Office  of  Environment  and  Energy;  800  Independence 
Avenue  S.W.;  Washington,  D.C. 20591 
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After  the  primary  register  has  been  loaded  the  secondary  register  is 
loaded  bv  depressing  the  following  kevs:  "t"  MP-SM*  300  STO  0  200  STO 
4  16  STO  6  l  CHS**  ATO  7  1.3  STO  8  "i”  "P-S* . 

Program  Exocut  ion 

The  program  is  started  bv  depressing  the  "E"  kev  for  the  assumed  "E" 
stability  (The  "B",  "C",  and  "D"  kevs  will  start  the  program  for  "Bn, 

"C",  and  "P11  stability  classes  respectively).  The  resulting  three 
numbers  printed  out  on  the  HP-97  or  displayed  on  the  HP-67  after  eaeh 
iteration  is  completed  are;  (1)  time  (in  seconds  )  from  the  program 
start;  (2)  distance  (in  meters)  that  the  point  in  the  emission  tail  has 
moved  down  the  runwav;  and  (3)  total  dose  (in  ppm-seconds )  that  the 
receptor  has  received. 

ft  is  noted  that  after  20  iterations,  the  dose  value  will  reach  a 
maximum  of  82.16  ppm- sec.  This  value  represents  the  dose  received  at  the 
receptor  from  the  first  emission  release  point  in  the  "tail".  When  the 
dose  converges  on  this  maximum  value  (when  all  concentration  digits 
remained  unchanged  out  to  the  second  decimal  point)  the  R/S  kev  is 
depressed  to  stop  the  program.  The  "A"  key  is  then  depressed  to  clear 
registers  and  index  the  starting  point  to  the  second  tail  position.  The 
"F."  key  is  then  depressed  a  second  time  to  start  the  next  computation. 

Again  when  the  dose  value  levels  off  at  55.59  ppm-sec.,  the  "R/S"  key  is 
depressed  to  stop  the  program.  Depressing  the  "A"  key  and  then,  after 
the  display  stops  flashing,  the  "E"  kev,  permits  the  last  computation  to 
be  completed  which  results  in  a  dosage  of  40.85  ppm-sec.  for  the  last 
tail  point.  After  reaching  this  last  convergence  value,  the  program  is 
terminated  bv  depressing  t he  R/S  kev.  Tl  '  person  making  the  calculation 
can  then  sum  the  three  dose  values  and  divide  them  by  3600  to  produce  a 
one-hour  concent ra t ion  of  0.05  ppm. 

CONCLUSIONS 

The  method,  limitations  and  use  of  the  SIMPLEX*^"  model  have  been  described. 
Fite  program  can  determine  concentrations  from  departing  aircraft  and  has 
tlu*  flexibility  to  easily  accept  parameter  changes.  Tt  can  treat  either 
single  or  multiple  events  and  permits  air  quality  calculations  to  be 
made  b\  persons  without  an  extensive  computer  background.  The  model  can 
assist  in  determining  the  impact  of  aircraft  emissions  on  air  quality  in 
conjunction  with  requirements  for  controlling  engine  emissions  and  can  be 
nst‘d  .is  a  screening  tool  in  evaluating  the  air  quality  impact  of  proposed 
Federal  actions  at  airports. 


*  The  P-S  command  is  input  by  depressing  the  "CLx"  key  on  the  HP-97  and  the 
"CHI”  key  on  the  HP-67  calculator. 

**  Negative  numbers  are  entered  into  the  Hewlett  Packard  calculators  by 
depressing  the  appropriate  number  key  followed  by  the  "CHS"  kev. 
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REGISTERS  AND  LABELS 


REG  I STERS 
Pr  imarv 

0  emission  rate 

1  'V  distance 

2  initial  sigma  ,,zn 

3  zero  register 

4  sigma  MzM  exponent 

3  sigma  My"  coefficient 

6  sigma  "y"  exponent 

7  sigma  "z"  coefficient 

8  total  distance  (1+3) 

9  sigma  "z" 

A  sigma  Myn 

B  plume  height 
C  "y"  distance 
D  wind  angle 
E  iteration  time 
I  wind  velocity 

Secondary 

0  fixed  source  receptor  distance  along  runway 

1  plume  rise  factor 

2  dose  summation 

3  sidewind  "y"  factor 

4  fixed  distance  -receptor  to  runway 

5  variable  distance  between  source  and  receptor  in  the  runway 

d  tree tion 

6  initial  sigma  "y" 

7  time  at  runway  location 

8  acceleration 

9  hypotenuse  ("h"  in  Figure  2) 

Ij\BELS 

A  program  to  clear  registers  and  make  required  inputs  for 
new  iteration  sequence 

B,C,D,E,  Storage  of  coefficients  and  exponents  for  sigma 
calculations  -  stability  classes  B,C,D,and  E 

4  subroutine  for  Label  A 

5  subroutine  for  Label  A 

6  subroutine  for  Label  A 

7  subroutine  to  switch  registers 

8  main  program  to  move  airplane  along  runway  and  to 

calculate  dosage 

9  subroutine  of  Label  8 


Figure  4 


INPUTS  and  OUTPUTS 


INPUTS 


Item  Units 

Primary  Register  (P) 


source  emissions  over 
duration  of  event- 
emission  rate  x  iteration  time 
initial  sigma  "z" 
plume  height 
wind  angle 
iteration  interval 
wind  velocity 


grams 

meters 

meters 

degrees 

seconds 

meters  per  second 


Secondary  Register 

fixed  source  receptor  distance 
along  runway 

fixed  distance  from  receptor 
to  runway 
initial  sigma  "y" 
beginning  time 
acceleration 


meters 

meters 

meters 

seconds 

meters/ sec/ sec 


OUTPUTS 


total  elapsed  time  at 
iteration 

fixed  source-receptor  distance 
along  runway 
dose  sum 


seconds 

meters 

parts  per  million  sec. 


Keys 


Sto  0 
Sto  2 
Sto  B 
Sto  D 
Sto  E 
Sto  I 


Sto  0 

Sto  4 
Sto  6 
Sto  7 
Sto  8 


Prnt  7 

Prnt  0 
Prnt  2 
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EMISSION 


TAIL  geometry 


Figure  6 


SAMPLE  PROBLEM  INPUTS  -  Case  1 


INPUTS 


1 1  em 


No./  Units 


Address  Primary  Register 

source  emissions  **over 
duration  of  event 
(emiss.  rate  x  iter,  time) 


Address  Secondary  Register 

fixed  source  receptor  distance 


Keys5* 


none 


(  80  gm/s  x  1  sec. 

) 

80 

grams 

Sto 

0 

initial  sigma  "z” 

8 

meters 

Sto 

2 

plume  height 

4 

meters 

Sto 

B 

wind  angle 

30 

degrees 

Sto 

0 

iteration  interval 

1 

second 

Sto 

E 

wind  velocity 

5 

meters  per  second 

S  to 

1 

f.  P-S 


along  runway 

+  300 

meters 

Sto 

0 

fixed  distance  from  receptor 
to  runwav 

200 

meters 

Sto 

4 

initial  si gma  "y" 

16 

meters 

Sto 

6 

beginning  time  *** 

-1 

second 

Sto 

7 

acceleration  **** 

1.3 

meters/ sec/ sec 

Sto 

8 

Readdress  Primary  Register 


*  Applicable  to  both  HP-97  and  HP-67  calculators  except  that  the 

wind  velocity  is  loaded  into  the  HP-67  calculator  by  depressing 
the  black  "hM  kev  followed  by  the  black  lettered  "ST  I”  key. 

**  Possible  data  source  -  (7). 

***  For  a  beginning  time  of  zero,  the  negative  value  of  the  iterative 
duration  must  be  input.  This  is  accomplished  by  entering  the 
duration  value  followed  by  the  "CHS"  key. 

Possible  data  source  -  Aircraft  Performance  Manuals. 


INPUTS 


Figure  7 

SAMPLE  PROBLEM  INPUTS  -  Case  2 

No./  Units 


Address  Primary  Register 

source  emissions  **over 
duration  of  event 
(emiss.  rate  x  iter,  time) 


Address  Secondary  Register 

fixed  source  receptor  distance 


Kevs* 


(  80  gm/s  x  1  sec.  ) 

80 

grams 

Sto 

0 

initial  sigma  "z" 

8 

meters 

S  to 

2 

plume  height 

4 

meters 

Sto 

B 

wind  angle 

30 

degrees 

Sto 

D 

iteration  interval 

1 

seconds 

Sto 

E 

wind  velocity 

5 

meters  per  second 

S  to 

I 

along  runway 

fixed  distance  from  receptor 

-300 

meters 

S  to 

0 

to  runway 

100 

meters 

Sto 

4 

initial  sigma  "y" 

16 

meters 

Sto 

6 

beginning  time  *** 

-1 

seconds 

Sto 

7 

acceleration  **** 

1.3 

meters/ sec/ sec 

Sto 

8 

Readdress  Primary  Register  f,  P- 

*  Applicable  to  both  HP-97  and  HP-67  calculators  except  that  the 

wind  velocity  is  loaded  into  the  HP-67  calculator  by  depressing 
the  black  "h"  key  followed  by  the  black  lettered  "ST  I"  key. 

**  Possible  data  source  -  (7). 

***  For  a  beginning  time  of  zero,  the  negative  value  of  the  iterative 
duration  must  be  input.  This  is  accomplished  by  entering  the 
duration  value  followed  by  the  "CHS"  key. 

****  Possible  data  source  -  Aircraft  Performance  Manuals. 


receptor 


wind  -  5  m/s 


runway 


airplane 

source 


negative 

(-262.5) 


Figure  8  (CONT.) 
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RESULTS  -  Case  2 
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figure  9  (Co NT.) 


.J  :• .  0  l* 

-275.  46 

ut 

r *  tt 

-2:5.35 

1*4 

^  ■.*  *  t  i 

44* 

al  1  ■  t  ; 

+  ♦1 

w  .  .  tl 

44* 

-  '  '?  7- 

III 

174.75 

444 

w  •  •  C'  1 

444 

_  J .  v  t- 

t*4 

-175.46 

m 

72.  60 

444 

215.  76 

444 

27.00 

m 

—  .  r 

_  I  •  TV 

444 

-127.  :z 

m 

77.  tiu 

441 

:-\ee 

i*i 

ir.B  5 

444 

_  -»r  JT. 

«  w  t 

Hi 

.  .  t 

444 

25.66 

Hi 

7^.66 

444 

-43. 75 

Hi 

??».40 

444 

2. 621*72754-65 

444 

v  <“.* .  C  L 

444 

27.  CO 

♦  H 

75.66 

4*4 

-10. €0 

i»» 

347.25 

444 

2.022 51555:6: 

HI 

7-:.  3l 

44* 

:7.w 

*H 

77. 66 

444 

-  7  or 

HI 

352. 46 

444 

-  'r?i'  ;  ?cr  *  _  -  • 

,  .  ICwCCw..  •.  -! 

444 

:  . .  w  i 

444 

Hi 

55.  7C 

Hi 

;.CT 

Hi 

25.66 

4*1 

55 

4*1 

*  .*  •* 

Hi 

75.0fi  444 

135.  444 

5.  rj  in 


REFERENCES 


1.  Wang,  I.T.,  L.  A.  Conley,  0.  M.  Rote,  Airport  Vicinity  Air 
Pollution  Model  Users  Guide,  Argonne  National  Laboratory, 
Argonne,  IL,  FAA  Report  No.  FAA-RD-7 5-230 ,  Dec.  1  975. 

2.  Slade,  D.  H.,  Meteorology  and  Atomic  Energy  1968;  U.S.  Atomic 
Energy  Commission,  July  1968. 

3.  Turner,  D.  8.,  Workbook  of  Atmospheric  Dispersion  Estimates, 

U.S.  Department  of  Health,  E  d  u cat  ion  a nd  We  1  f  a  r  e , 

Cincinnati,  Ohio,  Revised  1970. 

4.  Smith,  D.G.,  R.J.  Yarmartino,  C.  Benkley,  R.  Isaacs,  J.  Lee, 

D.  Chang,  Concorde  Air  Quality  Monitoring  and  Analysis  Program 
at  Dulles  International  Airport,  Federal  Aviation 
Administration,  Washington,  D.C.,  Report  No.  FAA-AEQ-7  7- 14 , 
December  1977. 

5.  Yarmartino,  R.J.,  D.G.  Smith,  S.A.  Bremer,  D.  Heinold, 

D.  Lamech,  and  B.  Taylor,  Impact  of  Aircraft  Emissions  on  Air 
Quality  in  the  Vicinity  of  Airports,  Federal  Aviation 
Administration,  Washington,  D.C.,  Report  No.  FAA-EE-80-09A&09B , 
July  1980. 

6.  Tank,  W.C.,  B.K.  Hodder,  '‘Engine  Exhaust  Plume  Growth  in  the 
Airport  Environment,  Proceedings  Air  Quality  and  Aviation:  An 
International  Conference,  Reston,  Virginia,  October  16-18,  1978, 
Federal  Aviation  Administration,  Washington,  D.C., 

Report  No.  FAA-EE-78-26 . 

7.  Segal,  H.W. ,  Emissions  from  Queuing  Aircraft,  Air  Pollution 
Control  Association,  73rd  Annual  Meeting,  June  1980, 

Montreal,  Canada,  Paper  80-3.5 

8.  U.S.  Environmental  Protection  Agency,  Compilation  of  Air 
Pollution  Emission  Factors  -  Third  Edition,  AP-42  Supplement  10, 
Research  Triangle  Park,  North  Carolina  27711,  February  1980. 


